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The Frenkel Kontorova (FK) model is known to exhibit the so called Aubry's transition 
which is a jamming or frictional transition at zero temperature. Recently we found similar 
transition at zero and finite temperatures in a super-conducting Josephson junction array 
(JJA) on a square lattice under external magnetic field. In the present paper we discuss how 
these problems are related. 

§1. Introduction 

Understanding of non-crystalline solids such as glasses and granular systems is 
an important problem in condensed matter physics)^'' 2 -' A useful concept is that 
frustration of geometrical, energetic or kinetic origins is indispensable to avoid crys- 
tallization and allow realization of amorphous solids.®'® In the present paper we 
discuss a jamming in a strongly frustrated Josephson junction array (JJA) under 
external magnetic It is a very interesting system which provides an excep- 

tional opportunity to study both athermal (jamming) and thermal (glass) transitions 
in exactly the same settings. The question raised by the Chicago group - whether 
athermal and thermal jamming or glass transitions can be understood in a unified 
- can be asked explicitly in this system. 

In the present paper we discuss the possibility that athermal and thermal jam- 
ming transition in the present system can be understood as a generalization of the 
Aubry's transitiorP^'^ found in a family of one-dimensional models of frictions, 
most importantly the Frenkel-Kontorova (FK) model which exhibits very rich phe- 
nomenology in spite of its simplicity!-^ 1 

The organization of the paper is as follows. In the next section, we discuss 
the sequence of connections between the FK modelp® Matsukawa-Fukuyama (MF) 
modeP® and the frustrated Josephson junction array under magnetic 

nekPHEJ 

step by step. In sec. we review the Aubry's transition^ 1 in the FK model. 
There we focus on the properties of the so called hull function which is a powerful 
theoretical tool to analyze the Aubry's transition. Then we sketch our recent attempt 
to generalize it for the case of frustrated JJA.® In seed] we point out that 'shear' can 
be exerted on JJA via external electric current.^ 1 We discuss how tribology (sliding 
friction),^ non-linear rheology (soft-matters, granular matters, e.t.e.}^H25l and 
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non-linear transport (JJA, superconductors, e.t.cp^*^* are related to each other 
emphasizing remarkable similarity of their scaling features around critical points 
including the J (Jamming)-point. Finally we discuss the "Jamming phase diagram" 
of the JJA, which is analogous to the one proposed for soft-matters pQJ3 suggested 
by our analysis of non-linear transport properties at zero temperature® and Monte 
Carlo simulations at finite temperatures.® In sec. [5] we summarize this paper and 
discuss some future outlooks. 

§2. Link between the friction models and the Josephson junction arrays 



Fig. 1. (Color online) Schematic pictures of the friction models and the Josephson junction ar- 
ray, a) the Frenkel-Kontorova (FK) model b) Matsukawa-Fukuyama's 2-chain model on the 
ladder lattice and c) the Josephson junction array (JJA) on a square lattice under external 
magnetic field. The squares in the plaquette represent positions of the vortexes induced either 
by mismatching between the two incommensurate surfaces b) or external magnetic field c). 



2.1. Frenkel-Kontorova model - starting point 
2.1.1. Frustration due to mismatching 

The original Frenkel-Kontorova model 12 ' is a one-dimensional elastic chain of 
particles put on a periodic substrate (see Fig. Q]a)). The Hamiltonian is given by, 



H 



n=l ^ 



U r , 



I) -Acos(— ) 



(2-1) 



Here u n denotes the position of the ra-th particle. The particles are connected to 
each other by Hookian springs of strength k as described by the 1st term in the 
Hamiltonian where I is the natural spacing between the particles. We impose a 
boundary condition such that the length of the whole system is fixed {ul + \ —u\)/l = 
The 2nd term describes the periodic potential due to the substrate whose period 



is a. 



*' The FK model with fixed volume (length) and that under fixed pressure (external force) 
behave completely differently EU>E2l q^hg f ormer j s relevant in the context of friction (jamming) and 
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There are two important parameters: 1) the strength of the potential A and 2) 
winding number /, 

' - b < 2 ' 2 > 

Both are crucial for the jamming-unjamming transition (Aubry's transition) in the 
FK modelPJ-TO Later we will find equivalent two parameters in the frustrated 
Josephson junction array under external magnetic field. 

The sinusoidal potential allows the elastic chain to make phase slips with respect 
to the substrate. Thus this simple system allows both elastic and plastic deforma- 
tions. The elastic term prefers to keep the natural spacing I while the substrate 
potential prefers 2a. In the context of the friction between two different materials 
brought in contact with each otherpS it is natural to suppose that the two surfaces 
are incommensurate with respect to each other, namely / is an irrational number - 
a number which cannot be represented as a ratio of some two integers. As the result 
the system becomes frustrated as soon as A becomes finite. Finding the ground state 
of the system, which is some periodic (possibly of very long periodicity) crystalline 
structure for rational /, becomes a hi ghl y non-trivial problemP^J In the present 
paper we always assume / is irrational. l2| 

The system exhibits a jamming or frictional transition - called as Aubry's tran- 
sitioiP™ which we review in section [3l For weak enough coupling A < A c the 
elastic chain is only mildly deformed and it can slide over the substrate smoothly 
without energy dissipation - friction-less. For stronger coupling A > A c , the elastic 
chain becomes pinned by the substrate and friction emerges. We will find later that 
the parameter A, which plays a key role in the FK model, is equivalent to strength 
of anisotropy of the Josephson coupling in the Josephson junction array. 

2.1.2. Phase representation 

It is convenient to introduce a dimension-less "phase" variable n defined by 
u n = ^ [9 n + 27r(/ — l)n] by which the (dimension-less) Hamiltonian can be rewritten 
as. 

H = Yl {V 6n+1 2lT ^ ~ Xcos ( 9 n + 2vr/n)| (2-3) 

n=l ^ > 

where A is also made dimension-less by an appropriate rescaling. The boundary 
condition is such that (6l+i — ^i)/2vr = L is fixed. 

frustration is in some sense stronger than the latter. In the latter case the response of the system 
with respect to the increments of the external force exhibits devil's stair case singularities. Note 
also that the Hookian spring force, which arises due to the harmonic potential, does not endure the 
natural spacing I by itself. 

*' For technical reasons we wish to use periodic boundary conditions which cannot be compatible 
with irrational f. Thus in practice we use rational numbers which approximates a target irrational 
number. For instance we can take a series of integers p n with n = 1, 2, . . . from the Fibonacci series 
and construct a series of rational numbers p n -i/Pn which converges to / = (3 - \/5)/20.38196601... 
in the limit n — > oo. We consider systems with linear size L — p n so that we arrive at the target 
irrational number in the thermodynamic limit L — > oo. Note that it is easy to construct similar 
Fibonacci-like series for any irrational numbers which are solutions of some quadratic equations.^ 
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2.2. Matsukawa-Fukuyama model - a crucial intermediate step 
2.2.1. Phase model on ladder-lattice 

Matsukawa and Fukuyama considered a two-chain model in the context of fric- 
tion.^ Their idea is to allow the "substrate" in the FK model to deform elastically 
as well, which is certainly more realistic than the FK model in the context of tri- 
bologyP^ Following their idea, let us modify the FK model Eq. (|2-3p and develop a 
phase model defined on a two- leg ladder lattice shown in Fig(T]b). The hamiltonian 
is given by, 



H 



-A Y, <x*(0i-6j-torfni) (2-4) 



To simplify notations we relabeled the sites as i = 1,2, . . . , N whose position in the 
real space is given by (rtj, mj). In the two-chain model the index for the column takes 
values n = 1,2, ... ,L while that for the row (or layer) takes just two values m = 1,2. 
The sums are took over nearest neighbour pairs connected by displacement vector 
= (rii — nj,mi — rrij) which is either equal to (1,0) or (0, 1). 

2.2.2. Gauge invariance 

An important property of the system is gauge invariance which we explain below. 
Let us rewrite the hamiltonian Eq. (|2-4p as, 



H = E E cos ^) ( 2 - 5 ) 

ey=(l,0) ey=(0,l) 

with the phase difference 

0i 0j A, r (2-6) 
Here Ay is an anti-symmetric matrix A^ = —Aji which satisfy, 

E A i,i = 2 *f- ( 2 - 7 ) 

plaquette 

The sum X^piaquette ^® ^ directed sum over "bonds" along each "plaquette" in the 
anti-clockwise manner. 

It is easy to see that the original representation Eq. (|2-4p respect the condition 
Eq. (|2-7|) . The crucial point is that the phase differences ipy are invariant under 
gauge transformations; 

0i -»• Bi + 56i (2-8) 

Aij -> A^ + 56i - 59 j (2-9) 



Thus the hamiltonian Eq. ()2-4p is gauge- invariant. The condition Eq. (|2-7p itself is 
also gauge-invariant. 

In addition to the gauge invariance, the hamiltonian Eq. (|2-5p is invariant under 
A -> -A with 6 n — > —0 n . Furthermore also / — > 1 -\- f does not change the problem. 
So we only need to consider < / < 1/2 in the following. 
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2.3. Frustrated Josephson-junction array (J J A) under magnetic field 

2.3.1. Frustration due to external magnetic field 

The final step is just to 1) increase the number of legs of the ladder to build 

a 2-dimensional square lattice and 2) replace the intra-layer elastic couplings by 

sinusoidal couplings ((see Fig. Q]c)). Then we obtain the Josephson junction array 

on a square lattice under external magnetic field applied perpendicularly to the 
arrayEJ 

H = - cos (^y) - A Yl cos (^i) ( 24 °) 

eij=(l,0) e i3 =(0,l 

Here Oi is identified as the phase of the superconducting order parameter of 
the i-th site (superconducting island) . The sinusoidal couplings represent Josephson 
coupling between the superconducting islands. Now the potential A is identified as 
the vector potential due to external magnetic field B z applied along the +z direction. 
The parameter / which appears in Eq. (|2-7p is the number density of quantized flux 
lines / = a 2 B z /(j)Q where a 2 ,B z and c/>o are the are of the plaquette, strength of the 
magnetic and flux quantum. 

Let us emphasize that the two important parameters in the FK model, namely 
1) the parameter A and 2) the winding number / are inherited down to the the 
JJA. To conclude we finally arrived a Josephson junction array on a square lattice 
with anisotropic coupling - with anisotropy A - under external magnetic field B z = 
(4>o/a> 2 )f - with irrational number density / of fluxes per plaquette. In short, let us 
call such a system as irrationally frustrated anisotropic JJA. 

Quite interestingly it is actually possible to construct anisotropic JUL in labora- 
tory. The strength of the Josephson coupling depends, for instance, on the thickness 
of the junctions. Saito and Osad£p3 created anisotropic JJA with various A by 
controlling the thickness of the junctions in the lithography process. 

It may sound rather strange to consider the anisotropy seriously since it usual 
plays only minor roles. Not surprisingly, previous studies of irrationally frustrated 
JJA considered only isotropic systems A = 10 As we discuss later, it turned out in 
our recent studied that A is actually relevant for irrational Quite remarkably 

the isotropic point A = 1 turned out to be a critical point at zero temperature 
corresponding to A c of the FK model where a jamming transition analogous to the 
Aubry's transition takes place. By symmetry it is obvious that we only need to 
consider the case A > 1. 

2.3.2. Vortex - analogue of dislocation 

We mentioned above that the parameter / can be regarded as number density 
of quantized flux lines per plaquette. As we explain below, this is because the vector 
potential A due to the magnetic field induces vortexes of the phases 6>j. The point 
is that vortexes are quantized objects like dislocations in crystals. 

*-* We note however that Denniston and Tang^J studied the frustrated JJA on the ladder-lattice 
(with m = 1,2) (See Fig. Q]b)) and consider variation of the inter-leg coupling A. Their system 
is almost the same as the 2-chain model by Matsukawa and Fukuyama but the elastic intra-chain 
coupling in Eq. (|2-5p is replaced by a sinusoidal coupling. They found the Aubry's transition also 
exist in the frustrated JJA on the ladder. 
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Here it is convenient to define "charges" of the vortexes as, 

q * = h s &ij)=Pi~f ^ = ...,-2,-1,0,1,2,... (2-11) 

plaquette 

where s(x) is a saw-tooth like periodic function with period 2tt and s(x) = x in 
the range — ir < x < vr. By definition, the charge q takes only discrete values of 
the form p% — f with some integer and offset — / as shown above. Physically the 
integer pi represents the number of quantized fluxes (each carrying a flux quantum 
4>q) threading the i-th plaquette. Note also that the charge q is gauge invariant. 

The usefulness of the charge becomes manifested in the so called coulomb-gas 
mapping (see Chap. 9 of !29p ) in which continuous, elastic deformations ( "spin- wave" ) 
are integrated out to find effective hamiltonian % of the vortexes. The resultant 
system is essentially equivalent to a lattice-gas of electrostatic charges interacting 
with each other by the repulsive coulomb interactions, 

U = nGifij)* + G(6) i (2' 12 ) 

with fij = (rii — rij,mi — rrij). 

The interaction potential G{r) is the (static) Green's function of elastic defor- 
mations (spin-wave). In 2-dimension, it scales as G(r) oc log(|f|) for r » 1, Note 
that the anisotropy A in Eq. (|2-10p is simply reflected in anisotropy in G(r) such 
that with it is stronger into y-direction f||(0, 1) compared to x-direction f||(l,0) by 
factor A(> 1). 

The value G(0)(> 0) can be interpreted as the core energy of the vortexes. 
Because of the core energy, states with higher values of the vortex charges generally 
have larger energies and can be neglected at low temperatures. Since we only need 
to consider < / < 1/2 as noted in sec 12. 2. 2l it is sufficient to consider two values 
of the charges q = —f, 1 — /. In addition we assume the charge neutrality ^ qi = 
holds, which can be enforced by applying the periodic boundary conditions. As the 
result we find that a fraction / of the plaquettes carries a vortex p = 1 (or q = 1 — /) 
and the other fraction 1 — / carries no vortex p = 0(orq = — /). In Fig. [H the 
boxes in the plaquette represent the vortexes (p = 1). 

2.3.3. Vortex patterns in equilibrium - vortex liquid, crystal and glass 

Let us sketch briefly possible patterns of vortexes in equilibrium states at low 
temperatures. For clarity we discuss three cases 1) / = 2) / is rational and 3) / is 
irrational. 

If / = 0, the ground state of the system is trivial: the phase becomes uniformly 
ordered 0i = constant for all sites i. In such a ground state the vortex is absent 
everywhere Pi = {qi = 0). It can be regarded as a crystalline state (or ferromagnetic 
state). At finite temperatures, pairs of vortex (p = 1) and anti-vortex (p = —1) will 
be created leading to melting of the crystalline state by proliferation of the vortexes 
(and anti- vortexes) at some critical temperature T c . In 2-dimension, it takes place 
in a special way named as Kosterlitz-Thouless transition!® 
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If / is rational, i. e. / = p/q with some integers p and q, the system will have 
a period vortex lattice^ 1 which is analogous to periodically ordered structure of 
dislocations in the so called Frank-Kasper phasePD For example with / = 1/2, the 
charges exhibit a checkerboard like order in which the sign of the charges alternates 
along x and y-axis as q = 1/2, —1/2, 1/2, —1/2, ... (or p = 1, 0, 1, 0, . . .). We also 
note that the 'half- vortexes' which appears in the case of / = 1/2 is identical to the 
so called chirality in frustrated magnets .ESJ^ESJ 

In bulk superconductors formation of the vortex lattice is well known. The latter 
is a triangular lattice called as Abrikosov lattice.^ 1 On the other hand, the vortex 
lattices in JJA are formed on top of the underlying square lattice so that it is a 
super-lattice. Thus the vortex lattices in JJA are usually pinned by the underlying 
lattice of the JJA while those in the bulk pure superconductors are free to move 
around unless some pinning centers are present .23 'S3 

Starting from the FK model we are naturally lead to consider irrational f. Ap- 
parently the system cannot develop simple periodic vortex lattices with irrational 
/ so that finding the ground state becomes a highly non-trivial problem. Indeed 
JJA with irrational / - irrationally frustrated JJA - has been regarded as a system 
which possibly exhibit a glassy phase since a seminal work by T. HalseyP^ This is a 
quite intriguing possibility since it means emergence of a glassy phase with frustra- 
tion but without quenched disorder - at variance with the conventional spin-glasses 
and vortex-glasses (superconductors with random pinning centers) which involve 
quenched disorder pSJ'ESJ Disorder may be somehow self-generated in this system. 
Indeed equilibrium relaxations of the irrationally frustrated JJA were similar to the 
primary relaxation observed in typical fragile supercooled liquids.^ 

§3. Low lying states and Aubry's transition 

3.1. Hull function of the FK model 

Now let us turn to review the Aubry's transition found in the FK modeP^^^'^ 
and related friction models including the MF modelP^ 1 '® A remarkable feature of 
the FK model is that mathematically rigorous analysis of the low lying states is 
possible based on the fact that configuration of the energy minima (and maxima) of 
the system satisfies a recursion relation which is identical to the so called standard 
map well known in dynamical systems. 

It is known rigorously that the ground state of the FK model can be expressed 

9 n = 2irn + G(fn + a) (3-1) 

where G is a periodic function with periodicity 1, i.e. G(x + 1) = G(x) for any x. The 
function G(x) is called 'hull function' and describes distortion of the configuration 
of the elastic chain due to the substrate potential. It is important to note that the 
entire region < x < 1 becomes equally populated in the thermodynamic limit 
L — > oo for irrational f. 

Quite remarkably the phase a is arbitrary, meaning that there is a manifold of 
ground states which have exactly the same energy. Moreover it is known rigorously 
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G(x) 


(meta)stable states 


A < Ac(/) 
A > Ac(/) 


analytic 

non-analytic, with infinitely many discontinuities 


only the ground state 
infinitely many metastable states^ 



Table I. Changes of the low lying states by the Aubry's transition 

that there is a "phase transition" for irrational /, called 'transition by breaking of 
analycity' (or Aubry's transition), at a critical strength of coupling A c (/). The basic 
feature of the transition is summarized in the table 13.11 In Fig. [2] we show the hull 
function of the FK model constructed from numerically generated ground states at 
various A (see Ref^ 1 for the method). 

From a physical point of view, a significant consequence of the Aubry's transi- 
tion is the "frictional transition" between the sliding phase and jamming (pinned) 
phased 1 ' Let us sketch the essence of the reasoning in the following. 

For A < A c , starting from a ground state, one can find a continuum of states 
with exactly the same energy by varying a. The point is that they are all related 
to each other by some continuous displacements of the particles in the real space. 
This comes from the fact that G(x) is analytic for A < A c . Thus no external force is 
needed to slide the whole system - friction-less or sliding. 

Existence of the sliding becomes trivial if the elastic chain itself is replaced by 
a rigid body. In such an extreme case of friction between two incommensurate rigid 
bodies, the forces between them oscillates in the space with an incommensurate 
period so that the net force becomes cancelled out. The non-trivial point is that 
similar cancellation of the forces still happens even if the chain is allowed to deform 
elastically as long as the coupling A is sufficiently small. 

For A > A c , discontinuous points appear in the hull function G(x). It means that 




0.4 0.6 
x=fn-int(fn) 

Fig. 2. (Color online) Hull function G(x) of the FK model. Here 9 n — 2irn of the ground states 
are plotted against the "folded coordinate" [n] — fn — int(/n) to elucidate the hull function. 
A c = 0.9716.. I 
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variation of a require discontinuous movements of the particles in the real space. 
Thus the ground states are no-more connected to each other by sliding: the system 
prepared in the ground state has to go over some higher energy states (thus energy 
barriers) to reach another ground state. Thus the system is jammed. Now some 
finite strength of external force greater than a certain frictional force (yield stress) 
/yield (A) oc (A — A c )^ must be applied to the system to let it move (de-pinning) ^ 

Firov et aP^'^ has been able to find a hierarchy of exponentially large number 
of low lying states on top of the ground state in the jammed phase A > A c . This is 
a very interesting observation from the view point of the physics of glasses. How- 
ever, unfortunately the FK model is an one-dimensional system so that the Aubry's 
transition disappears at finite temperatures. 

The frictional transition and emergence of discontinuity in the hull function 
has also been found in the Matsukawa-Fukuyama's 2-chain modelP^ Now it is very 
natural to expect that these features will be inherited down to our J JA on the square 
lattice under magnetic field. The main message that we find here is that we should 
vary the anisotropy A and see what happens in the low lying states. 

3.2. Low lying states of the anisotropic JJA 

Let us now turn to the anisotropic irrationally frustrated JJA with A > 1. 
Examples of the real space configurations of the vortexes in equilibrium at a low 
temperature are shown in Fig. The most prominent feature is the stripe pattern 
of the vortexes which are regularly stacked into y-direction (stronger coupling) and 
undulated along the x-direction (weaker coupling). 5 '" 71 The formation of the stripes 
is reasonable because the repulsive interactions between vortexes are anisotropic if 
A 7^ 1 as we noted in sec 12.3.21 

There are two important observations. First, the undulated stripe pattern is 
frozen in time, i. e. the ergordicity is broken. The pattern of the undulation cannot 
evolve dynamically by usual relaxational dynamics once such a structure is estab- 
lished. This is simply because the stripes are perfectly stacked into the y-direction in 
a belt. At a first sight, the stripe patterns may look similar to those found, for ex- 
ample, in liquid crystals. But they are very different because usual stripes fluctuate 
dynamically.^ 

Second, there is a family of low lying states with different patters of the trans- 
verse undulation as shown in Fig. Apparently the ground state should have no 
transverse undulation. Very interestingly the energies of the different patterns of the 
undulation shown in Fig. [3] are very close to each other suggesting a gap-less band 
of undulated states. Thus these undulated states are all relevant in the equilibrium 
ensemble. This is manifested in the structure factor of the vortexes which exhibit 
Bragg peaks into q y direction but a power law tail into q x direction.®'^ 

This is a very peculiar state of matter. Is this a glass? "No", in the sense that 
it has Bragg peaks which one would not expect for a glass. "Yes" , in the sense that 
there are many states with different patterns of undulation, which is a self-induced 
disorder, and they are separated by energy barriers. 

In a sense, the prediction by Halsey^ 1 - that superconducting glass (without 
quenched disorder) in the JJA with irrational f - is realized. However we must keep 
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Fig. 3. (Color online) Undulated vortexes in the anisotropic irrationally frustrated JJA. a) Exam- 
ples of real space patterns of vortexes in thermal equilibrium at a low temperature and b) in 
nearby energy energy minima. The thermalized configurations shown in a) are obtained by 
performing Monte Carlo (MC) simulations at T = 0.2 on a system with / = 21/55 and A = 1.5. 
The equilibration is extremely hard in this system so that we used a MC method which com- 
bines the Metropolis method, over-relaxation method and exchange MC method to endure 
equilibration.^ The filled squares represent plaquettes with vortexes with charge p = 1. The 
configurations in b) are obtained by minimizing the energy by simple energy descent algorithm 
starting from the thermalized configurations shown in a) . The energies of the energy minima of 
the configurations l)-3) are E = -5072.14311, -5072.38582, -5072.34445 respectively. 

in mind that here we are considering anisotropic JJA with A > 1 instead of the 
isotropic JJA A = 1 studied in most of the previous works. 

Now let us examine the low lying state more closer. In the analysis of the ground 
state of the FK model, the hull function Eq. (|3-ip played a central role as we noted 
before. Since the JJA can be regarded as a 2-dimensional version of the FK model, 
we are naturally led to look for similar one which may describe the low lying sates 
of the JJA in a compact way. 

Because of the gauge invariance, let us focus on the gauge-invariant phase differ- 
ences across the Josephson couplings ipij = Q% — 6j — Aij defined Eq. (j2-6j) where i and 
j are nearest neighbours across a Josephson coupling which may be either along x or 
y-axis. As we discuss later ipij is directly related to the Josephson current sm(ipij), 
which is the analogue of stress field in rheology. 

In Fig. [5] we display the phase differences ipij at various sites i = (n,m) plotted 
against "folded coordinates" [n] = fn — int(/n) and [m] = fm-mt(fm) which takes 
values limited in the range < [n] < 1 and < [m] < 1. The purpose of this plot 
is to elucidate the hull function analogously to the case of the FK model shown in 
Fig. [2j Quite remarkably the plots in the panels c) and d) strongly suggest there 
is indeed an analytic hull function of the folded coordinate along the direction of 
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stronger coupling. On the other hand, the panels a) and b) suggest there are no such 
analytic hull functions along the direction of weaker coupling. 

Recently we found it is possible to obtain the hull functions analytically by 
performing a 1/A expansion starting from A = oo limit.^ It turned out that the 
transverse undulation is encoded in the "phase differences" between different columns 
which one can see in the pane-11 c) and d). 

The existence (absence) of analytic hull functions along stronger (weaker) cou- 
plings immediately implies sliding (jamming) of the vortexes. Starting from an 
energy minimum, a family of different states with exactly the same energy can be 
obtained through the operation [m] — > [m+a] along the direction of stronger coupling 
with varying phase shift parameter a. This amount to a unidirectional motion of the 
undulated vortex stripes into the direction of stronger coupling without changing its 
pattern, i. e. sliding. In contrast, no such operation is possible along the direction of 
weaker coupling, i. e. jamming. In the next section we discuss how these properties 
are reflected in physical observables associated with shear. 

The above observation implies the symmetric system with A = 1, on which most 
of the previous works have been dedicated, is actually very special. As we discuss 
later, the critical point corresponding to the Aubry's transition point is actually 
A c = 1 in the JJA at zero temperature T = 0. 
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Fig. 4. (Color online) Spatial configuration of gauge invariant phases in an energy minimum. Phase 
differences Vi,*+e x an d ipi,i+e y across Josephson junctions parallel to x and y axes at various 
sites i = (n, m) are plotted against "folded coordinates" . Different symbols in the panels a) and 
b) represent the phase differences at n = 1, 2, . . . , L along arbitrary chosen three different "rows" 
of the JJA. In the panels c) and d), the phase differences at m — 1,2, ... ,L along arbitrary 
chosen three "columns" are shown. The system size is L = 55. 



§4. Response to shear - shear by external electric current 



Jamming is nothing but onset of rigidity which can be detected by response 
against shear. In general shear is induced into the system through the boundaries. 
In rheology one can consider to apply some constant external shear-stress cr ex t on 
boundaries of systems under study. Very interestingly this is equivalent to put 
external electric current I ex t into a Josephson junction array as we explain in sec 
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14.21 We will find that stress a and shear rate 7 in rheology correspond precisely to 
current I (or current density J) and voltage V (or electric field E) in the transport 
problem of the driven JJA. 

To study rheology (or transport) one can control either the external stress a or 
shear rate 7. Here we choose to control the strain 7 so that we can discuss static 
and dynamic response to shear in the same set up. 

We put shear across the Josephson junction array along, say y-axis, in the follow- 
ing manner. First we fix the configurations of the phase variables 61 on the bottom 
(m = 1) and top (m = L) layers. Second we slightly change the boundary such that 
a uniform displacement is imposed on the top layer (m = L) Oi — > 6% + (L — 1)7 
while the bottom layer (m = 1) is left in the same fixed configuration. This amount 
to induce a gradient of phase (16/ dy = 7 along the y-axis. Clearly 7 corresponds 
to shear-strain of the usual sense. As the result some internal stress a, which is 
super-current running across the Josephson junctions (see below), will be induced in 
the system. To study rheology, we drive the top wall with a constant speed so that 
the strain 7 increase with a constant shear-rate 7. Let us also remark that shear 
d9/dy = 7 on the phases along the y-axis amount to motion of vortexes (disloca- 
tions) into the orthogonal direction, i.e. x-axis. This is equivalent to say that the 
vortexes are driven by the Lorentz forced 

4.1. Static response to shear - static rigidity 

From static point of view, emergence of rigidity can be best quantified by shear- 
modulus. The free-energy of the system ^(7) can be formally expanded in power 
series of infinitesimal shear strain 7 as 



where a and fi are the shear-stress and shear-modulus respectively. (. . .} stands for 
a thermal average. 

Here 7 must be infinitesimal. The free-energy density F(j) /N, in the thermody- 
namic sense, must not depend on the boundary condition (including the shape of the 
container) so that shear-modulus must be zero in the thermodynamic sense even in 
solids. Thus when the shear-modulus [i defined by the fluctuation formula Eq. (|4-4p 
emerges, it means that the ordering of the 7—7-0 limit and the thermodynamic limit 
N — > 00 no more commute in sharp contrast to liquids. In turn this means linear 
elasticity must fail in solids. Physically this means that elasticity and plasticity must 
emerge simultaneously in solids.^ 

It is useful to note that the change on the boundary condition can be formally 
"absorbed" into the bulk part of the system by replacing the original Hamlitonian 



F{i) = F(0) + N{a)-f + -/i 7 2 + . . . 




Eq. dSITO]) by, 



H(r/) = - ^2 cos O;j)-A cos(ipij + 

gy=(l,0) eij=(0,l) 



(4-2) 
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Based on this observation we find that the stress a can be expressed as, 

Na= dmrd =x ^ gin( ^ (43) 

7 eij=(0,l) 

Similarly the shear- modulus // can be expressed as, 

H = b - [(a 2 ) - (a) 2 ] (4-4) 
where b is instantaneous or adiabatic shear-modulus ("Born term") defined as, 

m= &H(nl = x J2 cos(V^). (4-5) 

7 ey=(0,l) 

Fluctuation formulae for the elastic modulus like Eq. (j4-4p are well known in 
literaturePD In the context of XY models and super-conductors it is usually called 
as helicity modulus.^ The crucial term is the 2nd term which represents reduction 
of the shear- modulus due to thermal fluctuations of the stress a. 

In liquids, the two limits 7 — > and N — > 00 should commute. Then an identity 
/i = must hold meaning exact cancellation must take place between the Born term 
and the fluctuation term in Eq. (|4-4p . 

In the previous section we found the anisotropic irrationally frustrated JJA 
exhibits sliding/jamming in the low lying states such that the vortexes can slide 
freely along the stronger coupling but jammed along the weaker coupling. In turn 
this means that shear of the phases along stronger /weaker coupling causes finite/zero 
changes of the energy respectively. Consequently the shear-modulus fi must be 
finite/zero along stronger /weaker coupling at zero temperature T = 0. Indeed we 
observed this numerically.^ 

This is a quite intriguing situation - the anisotropic system A ^ 1 at zero tem- 
perature T = behaves either as solid or liquid depending on the axes along which 
one imposes the shear. From numerical observations it seems that the picture holds 
up to the symmetric point A = 1® suggesting that the symmetric point is actually 
the critical point A c = 1 where shear-modulus along a given axis becomes zero/finite. 

4.2. Dynamic response to shear - transport or rheology 

The shear-stress a defined in Eq. fj4-3|) is nothing but super-current flowing along 
y-direction in the Josephson junction array. More precisely according to the DC/ AC 
Josephson relations® the current Iij and voltage drop V{ — Vj across the junction 
are given by, 

Iij = sin(^) Vi -V 3 = ^L (4-6) 

Here we are assuming some appropriate rescalings to define the dimension-less quan- 
tities and Vi. 

At each site i (super-conducting island) the current must be conserved. By 
taking into account charging of the island and Ohmic energy dissipation we find, 

+ £ ^ir^ + £ ^ = «w - w) ( 4 -7) 

3 3 
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where the sums are took over nearest-neighbours. C and R are the capacitance of the 
islands and resistance of the junctions respectively. I ex t is the strength of external 
current which is injected from the top layer m = L and extracted from the bottom 
layer m = 1. Combining with the Josephson relation Eq. (|4-6p and the definition 
of the gauge- invariant phase difference ifiij given in Eq. (|2-6p . one easily finds an 
effective equation of motion of the phases 9i, which is called RCSJ (Resistively and 
Capacitively Shunted Junction) model. 141 Apparently it can be cast into the form 
of Newton's equation of motion, 

dOi dvi dH 

~dt =Vi m ^ + ^ + ^Z^\ V i- V i)= F ^\ 5 ra^L-d m ^l). (4-8) 

* 3 

which can be considered as a toy model for rheology of layered systems under external 
shear applied on the top and bottom wallsP^ 

From the above observations, it is clear that transport properties in JJA and 
rheology are quite analogous. Because of the shear, the velocity field d6i/dt will 
acquire a slope along the y-axis which can be identified with the shear rate 7. From 
the AC Josephson relation (the 2nd equation of Eq. (|4-6p ). we find that it mounts 
to a constant electric field E along the y axis. 

To summarize shear-stress a and shear-rate 7 in rheology correspond to electric 
current / (or current density J) and voltage drop V across the system (or electric 
field E) in the transport problem of JJA. Thus the so called "flow curves" in rheology 
corresponds to current-voltage IV (or JE) characteristics in JJA. In tribology we just 
need to consider only two layers m = 1,2 as in the Matsukawa-Fukuyama's 2-chain 
modelJ 1 -^ The yield stress o~ c is called as static frictional force. These problems have 
been studied extensively in the corresponding research communities but somehow 
the intimate analogy has not been appreciated.^ 1 

4.3. Non-linear rheology and transport 

Let us discuss here some basic phenomenological aspects of the non-linear rhe- 
ology and the non-linear transport associated with 2nd order phase transition, in- 
cluding the jamming transition. To be specific we will denote A — A c as the distance 
to the critical point which is natural in the context of the anisotropic JJA at zero 
temperature. However the readers can easily translate the discussion to different 
situations by replacing A — A c by distance to critical temperature T — T c or jamming 
density <j> — <j>j, e.t.c. depending on the problems at hand. 

Let us assume the following generic scaling form. 

° = i A " (ix^xp) < 4 ' 9) 

where /3 and A are critical exponents and the subscript ± stands for A > A c and 
A < A c respectively. Physically we expect the following behaviours: 1) Newtonian 
behaviour in the "sliding phase" (A < A c ), 2) Finite yield stress in the "jammed 
phase" (A > A c ) and 3) The explicit A dependence must disappear at the critical 
point (J-point) A = A c . Based on these intuitions let us conjecture the following 
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asymptotic behaviours of the scaling function cr(y), 



p±(y) 



y A < A c 

A>Ac (4-10) 
cy^ y > 1 



where <r(0) and c are some constants. Consequently the scaling ansatz predicts the 
following asymptotic behaviours (7 — > 0), 

f r?(A) 7 77(A) oc (A c - A)-(^-« A < A c 
lima = < c 7 ^ A = A c (4-11) 

I a c (A) cr c (A) = cr_(0)(A - Ac)^ A>A C 

Most importantly the power law fluid behavior a oc 7^/^ at the critical point is 
predicted. Usually f3/A < 1 which is called shear-thinning behaviour. 

For the transport problems in superconductors including JJA, one just need to 
replace shear-stress a by the electric current density J, shear-rate 7 by the electro- 
magnetic field E. The Newtonian law corresponds to the Ohmic law J = o~E with 
the linear conductivity and the yield stress a c corresponds to critical current J c . 

Recently the non-liner rheology of granular systems is found to obey this type of 
scaling around the J-pointPS^^D j n granular systems Bagnold's scaling must replace 
the Newtonian law in the unjammed phase. At least formally, the above argument 
can be easily modified to account for it. 

The above scaling ansatz is quite reminiscent of the scaling property of magne- 
tization of ferromagnetic models around the critical temperature T c . On purpose we 
actually used the same standard notations for the critical exponents, i. e. j3 and A, 
in the latter problem. Namely by replacing the stress a by magnetization m and 
strain rate 7 by magnetic field h, one recovers m oc \T — T c \^m±{h/\T — T C \ A ). One 
can easily find precise correspondences between 1) the Newtonian (Ohmic) law v.s. 
paramagnetic behaviour m = \h with the linear-susceptibility x diverging at T c 2) 
power law rheology a oc 7^/^ v.s. m oc h s with 5 = 13/ A at the critical points and 
3) yield stress (critical current) a c oc (A — A c )^ v.s. the spontaneous magnetization 
m s oc (T c - T)P. 

This type of scaling has been advocated first in the context of non-linear current- 
voltage characteristics of superconductors by Wolf, Gubser and IrnrypS They stud- 
ied non-linear current-voltage characteristics of superconducting film at the super- 
conducting phase transition, which is a Kosterlitz-Thouless type 2nd order phase 
transitionPS They pointed out the analogy with the scaling of the magnetization 
of ferromagnets. Such dynamical scaling ansatz has been extensively used in the 
studies of transport properties in high-T c superconductors, especially in the context 
of the vortex-glasses with quenched pinning centers.^ 1 

More recently Otsuki and SasgP^ has realized the same type of critical behavior 
in the context of the non-linear rheology of molecular glasses. Quite remarkably they 
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Fig. 5. (Color online) Flow curves of the FK model and JJA. Both models are simulated by the 
RCSJ method at zero temperature T = 0. The panel a) displays the master flow curve of the 
FK model using A c = 0.9716.. and /3 = A = 3. The panel b) displays the master flow curve of the 
irrationally frustrated anisotropic JJA using A e = 1, /3 = 1.19 and A — 3.5 so that /3/A = 0.34. 



were able to find a mean-field theory which predicts that the flow curves of the non- 
linear rheology are formally identical to the equation of state of the Landau-Ginzburg 
theory of ferromagnets under external magnetic field suggesting in particular j3/A = 
1/3- 

Let us now discuss the dynamical scaling properties of the FK model and the 
JJA under shear. We performed the RCSJ simulation on both models. The master 
flow curve of the FK model is displayed in the panel a) of Fig. [5] which follow well 
the expected dynamical scaling behaviour around the Aubry's transition point A c . 
A previous workP found a c oc (A - A c )^ with A c = 0.9716.. and 2.85 < ?jj < 3.06. In 
addition we found the system remains Newtonian for the entire sliding phase A < A c 
including the critical point so we assumed /3 = A in the scaling plot. 

For the irrationally frustrated anisotropic JJA, we pointed out in sec !4.1l that the 
shear- modulus [i along stronger /weaker coupling is finite/zero at zero temperature 
T = and that the symmetric point A c = 1 is the critical point where the shear- 
modulus [i along a given axes changes from finite/zero to zero/finite. Then it is 
quite natural to expect that the current-voltage curve of the system with respect 
to injection of the electric current along a given direction exhibit dynamical scaling 
feature at around A c = 1. This is indeed observed by a numerical simulation of 
the RCSJ dynamics.® The current- volt age curves collapse onto a master curves as 
shown in the panel b) of Fig. [5j 

4.4. Jamming phase diagram 

In Fig.[6]we show a schematic "jamming phase diagram" of the irrationally frus- 
trated anisotropic Josephson junction array. So far we discussed only the properties 
of the system at zero temperature T = and A > 1. Under the electric current J y 
injected along the y-axis, the system remains jammed as long as J y is smaller than 
the critical current J c ~ (A — . The configuration of the jammed solid phase is 
characterized by the frozen pattern of the undulated vortex stripes. Under strong 
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T 



T C (A) ~ (A - If 




vortex liquid 



vortex solid with 
frozen undulation 



J i, 



Jamming point 



Fig. 6. (Color onlie) " Jamming" phase diagram of the JJA 



enough current J y > J c , the system starts to move exhibiting the shear-thinning 
behaviour. On the other hand the same solid can slide freely with respect to electric 
current J x injected along the x-axis for the entire A > 1 region. For the region A <, 
we just need to interchange the x and y-axes in the above discussion. 

Previous studies on the irrationally frustrated JJA are almost exclusively con- 
cerned with the symmetric point A = 1. Recent intensive numerical studies at finite 
temperatures suggest T C (A = 1) = rJ33'H3 without finite temperature glass transi- 
tion anticipated in the early works On the other hand our recent studies on 
the static properties at low temperatures®'^ strongly suggest T C (A) > 1 at least 
sufficiently away from the symmetric point A = 1 and that T C (A) rapidly decreases 
approaching the symmetric point A = 1. These point toward the possibility of the 
jamming phase diagram depicted in FigO Quite interestingly it is very similar to 
the jamming phase diagram proposed by the Chicago groupP® Most notably the 
point (A = 1, T = 0) look quite similar to the Jamming point which deserves further 
studies. 

§5. Conclusions 

To conclude we discussed static and transport or rheological properties of the 
irrationally frustrated anisotropic Josephson junction array (JJA) which exhibits 
vortex stripes with self-generated randomness at low temperatures. We emphasized 
in particular the intimate connection between the friction models and the irrationally 
frustrated JJA which provides valuable insights into the problems. 



*' We only show the A > 1 part. Note that T c (l/A) = T C (A)/A holds due to the obvious 
symmetry between x and y axis. 
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